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Abstract 



In these talks, I discuss a few selected topics in integrable models that are of interest from 
various points of view. Some open questions are also described. 
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1 Introduction: 



The subject of integrable models now encompasses a very large area of research involving many 
seemingly different topics. It is not at all possible to give a detailed exposition of the subject in 
just a few lectures. Therefore, when the organizers of the Andre Swieca summer school asked me 
to choose a few topics on which to speak at the school, I agreed with a lot of trepidation. In order 
to make the lectures self-complete, I, of course, had to start with some basics. The subsequent 
topics that I talked about, naturally, represent a personal choice according to my interests. There 
are many other interesting areas that are being pursued by many active groups, but I simply could 
not have done justice to all, in the limited time available. Similarly, the literature on the subject 
is vast and it would have been totally impossible for me to even pretend to have a complete list 
of references. Consequently, I have only chosen a handful of references, for my talks, that I have 
absolutely used in the preparation of my lectures. I apologize to all whose works I have not been 
able to mention in my talks or references that I have not been able to list. 



2 Historical development: 



Let me begin with some introduction to the historical development of the subject. Let us begin 
with the simple, working definition of an integrable model as a physical system that is described by 
nonlinear partial differential equations, which can be exactly solved. We will make the definition 
more precise as we go along. There are quite a few systems of equations of this kind that arise in 
various physical theories. For example, in + 1 dimensions, the Toda lattice is described by the 
set of equations 

Qi = Pi, i = l,2,---,N 

= _e-('32-Qi) 

Pa = e~('^«~'5«-i) - e~('^«+i~'5«\ a = 2,S,--- ,N -1 (1) 

and consists of a chain of N particles on a one dimensional lattice at the coordinates Qi, with Pj 
representing the conjugate momenta. 



In 1 + 1 dimensions, similarly, there is the celebrated KdV equation described by 

where u, for example, may describe the height of a water wave from the normal surface. The 
variables of the KdV equation can be scaled to have arbitrary coefficients in front of all the terms. 
As a consequence, another form in which the KdV equation is also known corresponds to 

du du d^u 
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In 1 + 1 dimensions, there is also the non-hnear Schrodinger equation described by 



at dx^ 

+ 2k|V|'V^ (4) 



dt dx'^ 

where k is a constant measuring the strength of the nonhnear interaction and can be both positive or 
negative corresponding to an attractive or repulsive interaction. There are several other integrable 
systems in 1 + 1 dimensions, but these two are the ones that have been widely studied. There 
are fewer integrable systems in 2 + 1 dimensions, which include the Kadomtsev-Petviashvelli (KP) 
equations and the Davey-Stewartson (DS) equations. 



One of the most important features that all integrable models have is that they possess soliton 
solutions to the equations of motion. Solitons are defined as localized, non-dispersive solutions that 
maintain their shape even after being scattered. Historically, of course, research in this area grew 
out of J. Scott Russel's observation, in 1834, of a solitary wave travelling for miles maintaining 
its shape. It was only in 1895 that Kortewcg and dc Vries gave a mathematical description of 
such shallow water waves, which is known as the KdV equation. Being nonlinear and difficult to 
solve, these equations, however, did not generate a lot of interest. In 1965, Kruskal and Zabusky 
undertook a "computer" experiment, namely, they wanted to numerically study the evolution of the 
solutions of the KdV equation. What they found was impressive, namely, when certain solutions 
of the KdV equation were scattered off each other, they maintained their shape even after going 
through the scattering region. Kruskal coined the term "solitons" for such solutions in 1969 and 
the presence of such solutions generated an enormous interest in such systems from then on. 



Non-dispersive solutions: 



Most physical linear equations have dispersive solutions and the presence of non-dispersive 
solutions in such systems is quite interesting. To appreciate the origin of such solutions and to see 
their relation to the nonlinear interactions of the theory, let us analyze the KdV equation. 



du du d^u , , 

dt dx dx^ 



As we have noted, such a system describes shallow water waves, where we can think of u{x, t) as 
representing the height of the water wave from the normal surface of water. The first term on the 
right hand side represents the nonlinear term. Let us, for a moment, look at the KdV equation 
without the nonlinear term, namely, 

du d^u 

It is easy to write down the dispersion relation following from this equation, 

E{k) = (7) 
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This immediately tells us that the phase and the group velocities, associated with a wave packet, 
in this case, are different, namely, 



-^-e V -^-3k' (8) 

'-'phase — ^ — ) t^group — — \°J 

Thus, we see that if the KdV equation contained only the linear term on the right hand side, 
solutions will disperse. 

On the other hand, let us next assume that the KdV equation does not contain the linear term 
on the right hand side, namely, 

du du 

This is also known as the Riemann equation. This can be solved by the method of characteristics 
and the solution has the general form 

u{x,t) = f{x + ut) (10) 

which is quite interesting, for it says that the velocity of propagation is directly proportional to 
the height of the wave. Namely, the higher points of the wave will travel faster than those at a 
lower height. This is what leads to the breaking of waves etc. However, from our point of view, we 
see that this has a localizing effect, opposite of what the linear term leads to. The linear and the 
nonlinear terms, on the right hand side of the KdV equation, therefore, have opposing behavior 
and if they can balance each other exactly, then, we can have solutions that will travel without 
any dispersion. The presence of nonlinear interactions, therefore, is quite crucial to the existence 
of non-dispersive solutions. 

In the KdV equation, this indeed happens and we have non-dispersive solutions. For example, 
let us consider 



/v 

u{x,t) = Svsech^^ix + vt) (11) 



This gives 



— = -3v2sech'^^{x + vt)tanh.^{x + vt) 



du 5 A\Pv , \/v 

u— = — 9f2sech ——[x + vt)tanh——[x + vt) 

^ = _3^;fsech^^(a; + i;t)tanh^^(x + i;t) +6i;tsech^^(a; + vt)tanh^(x-|-i;t) 

With trigonometric identities, it is easy to check now that 

du du d^u 
dt dx dx^ 

holds so that this represents a solution of the KdV equation. From the form of the solution, it is 
clear that it maintains its shape as it travels (non-dispersive) and such a solution is known as a one 
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soliton solution. One can also construct multi-soliton solutions through what is known as Backliind 
transformations, which I will not go into. 



Conserved charges: 



Given the KdV equation, one can immediately construct three conserved charges, namely, it is 
easy to check that 




are conserved under the evolution of the KdV equation. Several people had also constructed up to 
13 conserved charges for the system when Kruskal conjectured that the KdV system has an infinite 
number of functionally independent conserved charges. This was subsequently proved and the 
conserved charges constructed through the Miura transformation (as well as through the method 
of inverse scattering). However, in retrospect, the presence of an infinite number of conserved 
charges associated with a system possessing soliton solutions is intuitively quite clear. As we have 
noted, solitons scatter through each other maintaining their shape. This implies that there must be 
conservation laws which prevents the solution from deformations. Since the soliton is an extended 
solution, there must, therefore, be an infinity of such conservation laws for the solution to maintain 
its shape through collisions. The important thing is that such a system has an infinite number of 
conserved charges, which also means that the system is integrable. 



Bi-Hamiltonian structure: 



Integrable systems are Hamiltonian systems. For example, in the case of the KdV equation, we 
note that if we define 

d 

{u{x),u{y)}i = d5{x -y) = -q^K^ - v) (13) 
then, the KdV equation can be written in the Hamiltonian form as 

du r , . „ -, d fl o d'^u\ du d^u ^, 
- = Mx), Hsh = ^ [-,u^ + 5^ j = + 5^ (14) 

However, what is even more interesting is the fact that the same set of equations can also be written 
in the Hamiltonian form if we define 

f 1/9 d \ \ 

so that 
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Namely, the KdV equation is Hamiltonian with respect to at least two distinct Hamiltonian struc- 
tures. 



Representing the two Hamiltonian structures as (operators acting on delta function) 

Vi = d 

^2 = d^ + \{du + ud) (17) 

we note that the first structure is what is normally called the Abelian current algebra, while the 
second structure is known as the Virasoro algebra. As a result, we do not have to worry about 
these structures satisfying the Jacobi identity, which they do. This is another general feature of 
integrable models, namely, the Hamiltonian structures of integrable models arc generally associated 
with symmetry algebras. In fact, some of the nonlinear algebras, such as the W algebras, were 
studied from the point of view of integrable models. 



Looking at the structure of the two Hamiltonian structures of the KdV equation, it is clear that 
not only does this system have two distinct Hamiltonian structures, but that that 

V = V2 + aVi (18) 

also represents a genuine Hamiltonian structure (not necessarily of the system) . This is a nontrivial 
statement, considering that a structure must satisfy the Jacobi identity - a nontrivial relation - in 
order to qualify as a Hamiltonian structure. In this case, it follows from the fact that T>2 is a 
genuine Hamiltonian structure for any variable u and that 

3 

V2{u + -a)=V2{u)+aVi (19) 

When two Hamiltonian structures have such a relation (namely, if two structures are Hamiltonian, 
then, a linear combination of the two is also), they are said to be compatible. When a Hamiltonian 
system can be described by two distinct Hamiltonian structures that are compatible, the system is 
said to be a bi-Hamiltonian system. 



The existence of two Hamiltonian descriptions for the same equation, of course, implies that 

du SH2 SH3 , , 

This is a prototype of the recursion relation that exists between conserved charges in such systems. 
One can define a recursion operator 

7^ = V^^V2 (21) 
which will relate the successive conserved charges as 

= 7^— ^ 22) 
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Furthermore, if the two Hamiltonian structures are compatible, one can further show that these 
conserved charges are also in involution with respect to either of the Hamiltonian structures, thereby 
proving that the system is integrable. 

The phase space geometry of integrable systems is quite interesting. These are, of course, 
symplectic manifolds, but because there are at least two distinct Hamiltonian structures (symplectic 
structures), these are very special symplectic manifolds. Let us call the two symplectic structures 
as ui and U2- Then, on this manifold, one can naturally define a nontrivial (1, 1) tensor as 

S = u^^uj2 (23) 

The evolution of this equation can be thought of as the Lax equation and, therefore, this gives 
a geometrical meaning to the Lax equation. Furthermore, one can show that if the Nijenhuis 
torsion tensor, associated with this (1. 1) tensor, vanishes, then, the conserved charges will be in 
involution. Consequently, the vanishing of the Nijenhuis torsion tensor can be thought of as a 
sufficient condition for integrability in this geometrical description. 



Let me also note here that since the KdV system has an infinite number of conserved quantities 
Hn, n= 1, 2, • • •, each of this can be thought of as a Hamiltonian and will lead to a flow as 

Thus, with every integrable system is a hierarchy of flows and these represent the higher order flows 
of the system. The entire hierarchy of flows shares the same infinite set of conserved quantities and 
are integrable. 



Initial value problem: 



An interesting question, in connection with these nonlinear integrable systems, is how can one 
solve the initial value problem. Namely, given the initial values of the dynamical variables, in 
such systems, how does one determine their values at any later time. In linear systems, we are 
familiar with techniques such as the Fourier transformation or the Laplace transformation, which 
help by transforming differential equations into algebraic ones. However, these methods are not 
very useful in dealing with nonlinear equations. The method that is useful (and, therefore, which 
can be thought of as the analog of the Fourier transformation in the case of nonlinear equations) 
is the method of inverse scattering. Let me explain this in some detail. 



Let us consider the linear Schrodinger equation 

(^d-' + ^uix,t)^^ = XiP (25) 

where d stands for ^ and u{x, t) is the dynamical variable of the KdV equation. Here t is just 
a parameter that the potential in the Schrodinger equation depends on and not the evolution 
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parameter of the Schrodinger equation. Since the potential u{x,t) depends on the parameter t, 
it follows that both ^ and A will depend on t as well. However, what Gardner, Greene, Kruskal 
and Miura observed was that, if u(x, t) satisfied the KdV equation, then, the eigenvalues. A, were 
independent of t, namely, the evolution is isospectral in such a case, or. 

At = (26) 

Furthermore, in such a case, the dependence of ^/^ on t is very simple, namely, 

V't = - {\u:, + a) V + f 4A + (27) 
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In such case, the evolution of the scattering data, such as the reflection coefflcient, the transmis- 
sion coefficient etc, with t is easy to determine and, in fact, take a simple form. Thus, the strategy 
for solving the initial value problem can be taken as follows. Let us choose the linear Schrodinger 
equation with the potential u{x, 0) and determine the scattering data. Next determine the scatter- 
ing data at an arbitrary value of t from the simple evolution of the scattering data. Once we have 
the scattering data for an arbitrary t, we can ask what is the potential, u{x, t), which would give rise 
to those scattering data. This is essentially the method of inverse scattering. The reconstruction of 
the potential from the scattering data is done through the Gel'fand-Levitan-Marchenko equation, 

roo 

K{x,y) + B{x,y)+ dz K{x, z)B{y + z) = 0, y>x (28) 

Jx 



1 POO ^ 

(x) = — y dk R{k)e'^'' + Cne"''"'' (29) 



where 

B 

9-7r / 

n=l 

Here, R is the coefflcient of reflection, «;„'s represent the eigenvalues for the bound states and c^'s 
correspond to the normalization constants for the bound state wave functions. Once the solution 
of the Gel'fand-Levitan-Marchenko equation is known, the potential is determined from 

1 dK{x,x) 

-u{x) = 2^— (30) 



Let us see explicitly how the method works, in an example. However, let me also note that, 

while one can, in principle, find a solution to the GLM equation, in practice it may be difficult 
unless the starting potential were very special. One such class of potentials are solitonic potentials 
which are known to be reflectionless. In such a case, 

= (31) 

and this makes calculations much simpler. 
Let us, therefore, consider 

u{x, 0) = 12sech2a; (32) 
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We recognize, from our earlier discussion, that this is the one soUton solution of the KdV equation. 
Such a potential leads to no reflection. It supports only one bound state, for which 

K= 1, tpixjO) = -sechx (33) 

Therefore, we obtain 

c(0) = (y^~ dxij'^{x,0)^ = 2 (34) 
Prom the equation for the "time" evolution of the wave function, it is easy to determine that 

c{t) = c(0)e-^* = 2e-^* (35) 

so that we have 

B{x, t) = c{t)e-'' = 2e-^*-^ (36) 
In this case, the GLM equation becomes, 

K{x, y, t) + 2e-^*-^-^ + 2 / dz K{x, z, t)e-^*-y-' = (37) 

Jx 

This determines 

2g-8t-a;-2/ 

K{x,y,t) = - 



1 + e- 



-8t-2x 



dKix x^ 

or, u{x,t) = 12 = 12sech2(a; + 4t) (38) 

ox 

This is, of course, the solution that wc had determined earlier (corresponding to a specific choice 
of v) and this explains how the method of inverse scattering works. (It is worth noting here that 
the inverse scattering method was also independently used by Faddeev and Zakharov to solve the 
KdV equation.) 

The Lax equation: 

In some sense, the Lax equation is a formal generalization of the ideas of Gardener, Greene, 
Kruskal and Miura. Let us consider a linear operator, L{t) that depends on a parameter t through 
the potential. Let us assume the eigenvalue equation 

L(t)V' = AV' (39) 

along with the evolution of the wave function 

(40) 

where B represents an anti-symmetric operator. It follows now that 

dL{t) , M ^ , Bilj 

dL 

or, —ip + LBip = Xt'ip + XBip = Xfip + BLip 
ot 

dL 

or, —V = [S,L]V + AtV (41) 
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It follows, therefore, that 

At = (42) 

provided 

§ = lB.il (43) 



This is known as the Lax equation and L, B are called the Lax pair. What this equation says is 
that the evolution of the linear equation with respect to the parameter t will be isospectral, provided 
the Lax equation is satisfied (A is commonly referred to as the spectral parameter.). Furthermore, 
if a Lax pair is found such that the Lax equation yields a given nonlinear equation, then, this says 
that one can associate a linear Schrodinger equation with it and the method of inverse scattering 
can be carried out for this system leading to the integrability of the nonlinear system. This is the 
power of the Lax equation and, as we have mentioned earlier, the geometrical meaning of the Lax 
equation is that it represents the evolution of a special (1, 1) tensor in the phase space (symplectic 
manifold) of the system. 

As an example, let us analyze the KdV equation in some detail. Let us note that if we choose, 



then, with some straight forward computation, we can determine that 



so that the Lax equation 



leads to 



L{t) = d'^ + -u{x,t) 
B{t) = Ad^ + ^{du + ud) 



Idu 1 I du d^u\ 

6 dt 6 [ dx dx^ I 

du du d^u 

' dt dx dx^ 



(44) 
(45) 



We recognize this to be the KdV equation and, having a Lax representation for the equation, then, 
immediately determines the linear Schrodinger equation associated with it, which we have described 
earlier, in connection with the method of inverse scattering. This is, however, a general procedure 
that applies to any integrable model and that is why the Lax equation plays an important role in 
the study of integrable systems. 
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Zero curvature formalism: 



There is an alternate method of representing integrable systems, which brings out some other 
properties associated with the system quite nicely. Let us continue with the example of the KdV 
equation and consider the following vector potentials. 

7 f Q^X + gCxx 18^^ "^5"^^ ^ 

- ( f ) (*^) 

There arc several things to note here. First, C = C[u,X\ and that the vector potentials belong to 
the Lie algebra of SL{2,K), namely, Ao,Ai G SL{2,K). 

The curvature (field strength) associated with these potentials can be easily calculated and 
one recognizes that the vanishing of the curvature yields the equations associated with the KdV 
hierarchy. This is seen as follows. 

Foi = doAi - dlo - [Ao, Ai] = (48) 

gives 

Ut = C^xx + lidu + ud) C - 4AC^ (49) 

For A = and C = u, this coincides with the KdV equation. In general, we can expand in a power 
series of the form 

N 



C = ^(4A)^-"C„M (50) 

n=0 

Substituting this into the equation and matching the corresponding powers of (4A), we obtain 

Co = 1 

(d^ + ^{du + ud)jc„ = dCn+i, n = 0,l,2,---,iV-l 

Ut = (d^ + ^{du + ud)^CN (51) 

We recognize these as giving the recursion relation between the conserved charges (which we have 

discussed earlier) as well as the th equation of the hierarchy. This is known as the zero curvature 
representation of the integrable system and brings out the recursion relation between the conserved 
charges, the current algebra etc quite nicely. 



Drinfeld-Sokolov formalism: 



Thus, we see that an integrable model can be represented as a scalar Lax equation as well as a 
matrix zero curvature condition. The natural question that arises is whether there is any connection 
between the two. 
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To analyze this question, let us note that the scalar Lax equation for KdV is described by the 
Lax pair which leads to the linear equations 

Li) = (^"^ + ^vL^i) = Xi; 

^ = = (^d^ + ^ (5u + ud)^ V (52) 

The scalar Lax equation can be thought of as the compatibility condition for these two equations 
when the spectral parameter is independent of t. We note that, while the second equation may 
appear to be third order in the derivatives, with the use of the Schrodinger equation, the higher 
order derivatives can, in fact, be reduced. As a result, this pair of equations appears to be at the 
most quadratic in the derivatives. 



We know that a second order equation can be written in terms of two first order equations. 
Keeping this in mind, let us define 

^1 = (5^) (53) 
as well as a two component column matrix wavefunction 

(«) 

It is clear now that the linear Schrodinger equation 



can be written in the matrix form as 

Similarly, the time evolution equation (depending on B) can also be written as a matrix equation 
of the form 

dt^ = Ao* (56) 
The compatibility of these two matrix linear equations leads to the zero curvature condition 

dtA^-d^Ao-[Ao,Ai] = Q (57) 

These potentials, however, do not resemble the potentials that we studied earlier. However, it is 
easy to check that the two sets of potentials are related by a global (A is a constant independent 
of t) similarity transformation. For example, note that 

A^ = S-'A^S, ^={ \ "I^) (^^) 

This, therefore, establishes the connection between the two formalisms and tells us how to go from 
one to the other or vice versa. 
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3 Pseudo-differential operators: 



With the basics of the previous section, we are now ready to discuss some of the integrable models 
in some detail. The first thing that we note is that the Lax formalism and the Lax pair is quite 
crucial in the study of integrable models. However, finding a Lax pair, for a given integrable model, 
seems like a formidable task. This is where the Gel'fand-Dikii formalism comes to rescue. 

Let us consider a general operator of the form 



If z > 0, namely, if the operator P only contains non-negative powers of d, then, it is a differential 
operator (i = term is a multiplicative operator) On the other hand, if i takes also negative values, 
then, the operator P is known as a pseudo-differential operator (Formally, is defined from 
dd~^ = 1 = d~^d). There are some standard nomenclature in using pseudo-differential operators. 
Thus, for example, 




(59) 



where 




Oj = ai{x) 



(60) 




(61) 



and, correspondingly. 




(62) 



By construction, therefore, we have 



P = P+ + P- 



(63) 



We can also define, in a corresponding manner. 




i>k 



(64) 



Let us now note the standard properties of the derivative operator, namely. 




(65) 



which holds true for any positive or negative, where 




(66) 
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and Z*^*^) denotes the k th derivative of the function /. It is worth noting here from the above 
formulae that positive powers of the derivative operator acting to the right cannot give rise to 
negative powers of the derivative and vice versa. 

Using these properties of the derivative operators, we note that we can define a multiphcation 
of pseudo-differential operators. The product of two pseudo-differential operators defines a pseudo- 
differential operator and they define an algebra. We can define a residue of a pseudo-differential 
operator to be the coefficient oi d~^, in analogy with the standard residue, namely. 

Residue P = Res P = a_i(x) (67) 

This allows us to define a concept called the trace of a pseudo-differential operator as 

Trace P = Tr P = j dxRes P = jdxa-i{x) (68) 

Let us note that, given two pseudo-differential operators, P, P', 

Res [P,P'] = {df{x)) (69) 

In other words, the residue of the commutator of any two arbitrary pseudo-differential operators is 
a total derivative. This, therefore, immediately leads to the fact that 

Tr PP' = Tr P'P (70) 

since the "trace" of the commutator would vanish with the usual assumptions on asymptotic fall 
off of variables. This shows that the "Trace" defined earlier satisfies the usual cyclicity properties 
and justifies the name. 

Let us also note that, given a pseudo-differential operator 

i 

we can define a dual operator as 

Q = Y,d-'Q-i (72) 

i 

where, the g^'s are independent of the a^'s. This allows us to define a linear functional of the form 

Fq{P) = IV PQ = j dxY, mi-i (73) 



The Lax operators, as we have seen earlier in the case of the KdV equation, have the form of 
differential operators. However, in general, they can be pseudo-differential operators. Thus, there 
exist two classes of Lax operators. Operators of the form 

L„ = a" + uia"-i+ii25"-2 + ... + «„ (74) 
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are differential operators and lead to a description of integrable models called the generalized KdV 
hierarchy. On the other hand, Lax operators of the form 



(75) 



correspond to pseudo-differential operators and lead to a description of integrable models, com- 
monly called the generalized KP hierarchy. 

Let us consider the Lax operator for the generalized KdV hierarchy, for the moment. Thus, 

Ln = d^ + + • • • + Un (76) 

We can now formally define the n th root of this operator as a general pseudo-differential operator 
of the form 

oo 

{Ln)'^ =d + Y,ai{x)d-' (77) 



such that 



j_\ n 
Ln 1 = Ln 



(78) 



This allows us to determine all the coefHcient functions, ai{x), iteratively and, therefore, the n th 
root of the Lax operator. 



Let us next note that, since 



for any k, it follows that 

dLn 
dtk 



Ln 1 , Ln 



k 

J^n 



1 Ln 



k ^ mn 



(79) 



(80) 



defines a consistent Lax equation. This can be seen as follows. First, ii k = mn, then. 



T ' 



(L^^l 



and, therefore, the commutator will vanish and we will not have a meaningful dynamical equation. 
For k ^ mn, we note, from the structure of the first commutator, that it will, in general, involve 
powers of the derivative of the forms 

QU+k—l gn+k—2 ^0 

On the other hand, the terms in the second commutator will, in general, have powers of the 
derivative of the forms 

g„-2gn-3Qn-4___ 

However, if the two expressions have to be equal, then, they can only have nontrivial powers of the 
derivative of the forms 

Qn-2gn-3___g0 
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This is precisely the structure of the Lax operator (except for the term with which says that 

the above equation represents a consistent Lax equation. This equation also will imply that 

which is why, often, this constant is set to zero (as is the case in, say, the KdV equation). This 
result is very interesting, for once we have a Lax operator, the other member of the pair can now 
be identified with 

Bk = (d\ (81) 



up to a multiplicative constant. Such a Lax representation of a dynamical system is known as the 
standard representation. 



Furthermore, we note that the Lax equation also implies that 

dtk 

It is straight forward to show, using this, that 

dtA^Ln = dt^dt^Lr^ (83) 

Namely, different flows commute. This is equivalent to saying that the diff'erent Hamiltonians 
corresponding to the different flows are in involution. Therefore, if we have the right number of 
conserved charges, the system is integrable. 



T " 



T " 
, J-in. 



+ 



(82) 



The construction of the conserved charges, therefore, is crucial in this approach. However, we 
note from 



that 



dtk 



a 1 

—TV Ln=Ti 
dtk 



T ' 



T ' 



k 

J-in 



r n 



= 



(84) 



(85) 



which follows from the cyclicity of the "trace" . Therefore, we can identify the conserved quantities 
of the system (up to multiplicative factors) with 



m 



T n 
J-in 



m ^ In 



(86) 



This naturally gives the infinite number of conserved charges of the system, which, as we have shown 
before, are in involution. Therefore, in this description, integrability is more or less automatic. 
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Hamiltonian structures: 



In the Lax formalism, we can also determine the Hamiltonian structures of the system in a 
natural manner. These are known in the subject as the Gel'fand-Dikii brackets and they are deter- 
mined from the observation that the Lax equation looks very much like Hamilton's equation, with 

^L^^ playing the role of the Hamiltonian and the commutator substituting for the Hamiltonian 

structure. Analyzing this further, one ends up with two definitions of Gel'fand-Dikii brackets, 
which give rise to the two Hamiltonian structures of the system. With the notation of the linear 
functional defined earlier, they can be written as 

{FQ{Ln),Fv{Ln)}i = Tr (L„[F,Q]) 

{FQ{Ln),Fv{Ln)}2 = TV (^LnQ{LnV)^-QLn{VLn)+) (87) 

The second bracket is particularly tricky if the Lax operator has a constrained structure and the 
modifications, in such a case, are well known and I will not get into that. 



It is worth noting that these brackets are, by definition, anti-symmetric as a Hamiltonian struc- 
ture should be. While the first bracket is manifestly anti-symmetric, the second is not. However, 
it is easy to see that the second is also anti-symmetric in the following way. 



{FQ{Ln),FviLn)} = TV [LnQ{LnV)^-QLniVLn) 

= Tr (LnQLnV -LnQ{LnV)_-QLnVLn + QLn{VLn)_) 

= Tr {-{LnQ)+{-LnV)_ + {QLn)+{VLn)_) 

= Tr (-{LnQ)+LnV + {QL^)^VLn) 

= -Tr (LnV{LnQ)^-VL^{QL^)^) 

= -{Fv{Ln),FQ{Ln)} (88) 

Thus, the two brackets indeed satisfy the necessary ant-symmetry property of Hamiltonian struc- 
tures. Furthermore, it can also be shown (I will not go into the details) that these brackets satisfy 
Jacobi identity as well and, therefore, constitute two Hamiltonian structures of the system. 



Without going into details, I would like to make some general remarks about the Lax operators 
of the KP type. Let us consider a Lax operator of the type 



(89) 



This is a pseudo-differential operator, unlike the earlier case, and, as we have already remarked, 
such Lax operators describe generalized KP hierarchies. In this case, it can be shown that a Lax 
equation of the form 



k 

An 
n 



>m 



(90) 
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is consistent, only for m = 0, 1,2. This is, therefore, different from the generahzed KdV hierarchy 
that we have aheady studied. For m = 0, the Lax equation is called, as before, a standard 
representation, while for m = 1, 2, it is known as a non-standard representation. All the ideas that 
we had developed for the standard representation go through for the non-standard representation 
as well and we will return to such an example later. 



Example: 



As an application of these ideas, let us analyze some of the integrable models from this point 
of view. First, let us consider the KdV hierarchy. In this case, we have already seen that 

L = L2=d^ + ^u (91) 

As we had noted earlier, we note that the coefficient of the linear power of d has been set to zero 
(which is consistent with the Lax equation). In this case, we can determine the square root of the 
Lax operator, following the method described earlier, and it has the form, 

=d + —ud-^ - —u^d-^ + —(u^x- -uA d-^ + --- (92) 
12 24 48 V 6 y ^ ^ 

It is easy to check that the square of this operator leads to L up to the particular order of terms. 



In this case, we have 
which gives 



1 

L2 



d 



(93) 



or. 



dL_ 
dti 
du 



L 



[d,L] 



du 
dx 



(94) 



This is the chiral boson equation and is known to be the lowest order equation of the KdV hierarchy. 
Let us also note that 



(lI) =(LLi)_ 



+ -ud + \u:c = d^ + \{du + ud) = \b 
4 o 8 4 



(95) 



where B is the second member of the Lax pair for the KdV equation that we had talked about 
earlier. It is clear, therefore, that 

dL 



3 

L2 



(96) 



will lead to the KdV equation. Similarly, one can derive the higher order equations of the KdV 
hierarchy from the higher fractional powers of the Lax operator. 



We note from the structure of the square root of L that 

1 

12 



Tr = — I dxu{x) (97) 
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Similarly, 



Tr lI 



1 2\ 1 2 

Uxx - 7:u ] + —u 



1 

96 



72 

dxu\x) (98) 



Up to multiplicative constants, these are the first two conserved quantities of the KdV hierarchy 
and the higher order ones can be obtained similarly from the "trace" of higher fractional powers of 
the Lax operator. 

Prom the form of the Lax operator, in this case, 

L = d^ + \u 
6 

we note that we can define the dual operators 

Q = d-^q2 + d-^qi, V = 8-^2 + d'^vi (99) 

Here gj and Vi are supposed to be independent of the dynamical variable u, so that the linear 
functionals take the forms 

Fq{L) = Tr LQ = ^J dxuqi, Fy (L) = ^ J dx uvi (100) 
In this case, we can work out 

{FQiL),Fv{L)}i = ^Jdx dyqi{x)viiy){u{x),u{y)h (101) 

On the other hand, 

TV L[V,Q]= J dx {qivi,^ - qi,^vi) = -2 J dxqi^^vi (102) 

Thus, comparing the two expressions, we obtain 

{u{x),u{y)}i = 72^5{x-y) (103) 

We recognize this to be the correct first Hamiltonian structure for the KdV equation (except for a 
multiplicative factor). The derivation of the second Hamiltonian structure is slightly more involved 
since the structure of the KdV Lax operator has a constrained structure (the linear power of d is 
missing). However, the construction through the Gel'fand-Dikii brackets, keeping this in mind, can 
be carried through and gives the correct second Hamiltonian structure for the theory. 

Let me note in closing this section that the generalization of the method of inverse scattering 
as well as the generalization of the Lax formalism (or the Gel'fand-Dikii formalism) to higher 
dimensions is not as well understood and remain open questions. 
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4 Two boson hierarchy: 



In this section, I will describe another integrable system in 1+1 dimensions, which is very interesting. 
The study of this system is of fundamental importance, since this system can reduce to many others 
under appropriate limit /reduction. It is described in terms of two dynamical variables and has the 
form 



du 

= {2uh + ah^)^ (104) 



— + « -auxj 
ot \ y X 



= (2h ' -2 

dh 
'di 



Here a is an arbitrary constant parameter and we can think of h as describing the height of a 
water wave from the surface, while u describes the horizontal velocity of the wave. This equation, 
therefore, describes general shallow water waves. 



This system of equations is integrable and, as we have mentioned, reduces to many other 
integrable systems under appropriate limit/reduction. To name a few, let us note that, when the 
parameter a = 0, this system reduces to Benney's equations, which also represents the standard, 
dispersionless long water wave equation. For a = —1 and /i = 0, this gives us the Burger's equation. 
When a = 1 and we identify 

u= , 11 = QQ 

q 

we obtain the non-linear Schrodinger equation from the two boson equations. Similarly, both the 
KdV and the mKdV equations are contained in this system as higher order flows. Thus, the study 
of this system is interesting because once we understand this system,, properties of all these systems 
are also known. 



Conventionally, the two boson equation is represented with the identifications 

a = 1, u = Jo, h = Ji (105) 

which is the notation we will follow in the subsequent discussions. In these notations, therefore, 
the two boson equations take the form 

= (2Jo Ji + Ji,a.)^ (106) 



Being integrable, this system of equations can be described by a Lax equation. Let us consider 
the Lax operator 

L = d- Jo + d-^Ji (107) 
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so that it is a pseudo-differential operator. Let us note that, for this operator, 

(l^)^^ = 53 - 3Joa2 + 3( Jl + 4 - Jo,x)5 
and so on. It is now straight forward to compute and show that 

L, = - (2J1 + Jl - Jo,x)^ + d-^ (2JoJi + Ji,^) 



It is, therefore, clear that the non-standard Lax equation 

'dt 

gives as consistent equations the two boson equations. 



L 



(108) 



(109) 



(110) 



The higher order flows of the two boson (TB) hierarchy can be obtained from 

dL_ 

dtk 

Of particular interest to us is the next higher order equation coming from 

dL_ 



L 



These have the forms 

dJo 
dt 
dJi 

'dt 



-Jo,xxx — 3 {JoJo,x)x ~ 6 {JoJl)x ~ ("'0)^ 

-{jl,xx + ^{JoJl)x + ^{jl{Jl + Jq — Jo,x 



(111) 



(112) 



(113) 



This equation is interesting, for we note that if we set Jq = and identify Ji = ^u, the equations 
reduce to the KdV equation. That is, as we had mentioned earlier, the KdV equation is contained in 
the higher flows of the TB hierarchy. Furthermore, it is also interesting to note that this provides a 
nonstandard representation of the KdV equation, unlike the earlier example where it was described 
by a standard Lax equation. 



Given the Lax representation, the conserved charges are easily constructed by the standard 
procedure from 

(114) 



so that we have 



Hn = Tr U 



Hi 


= Tr 




dx Jl 


H2 


= Tr 




j dx JqJi 


H3 


= Tr 


L3 = 


J dx (jf 



(115) 
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and so on. It is clear that if we set Jq = and identify Ji = all the even conserved charges 
vanish and the odd ones coincide with the conserved charges of the KdV equation. 



Hamiltonian structures: 



Let us denote the generic Hamiltonian structure associated with this system as 

(as {ti:})-"*'^-"' <^^^> 

Then, it can be shown that the TB equation has three Hamiltonian structures. But, let me only 
point out the first two here. 

so that we can write the two boson equations as 

It is now easily checked that, under Jq — ^ Jq + a, where a is an arbitrary constant, 

V2^V2 + aVi (119) 

which proves that these two Hamiltonian structures are compatible and that the system is integrable 
(which we already know from the Lax description of the system). 



Normally, the second Hamiltonian structure of an integrable system is related to some symmetry 
algebra. To see this connection, in this case, let us redefine (this is also known as changing the 
basis) 

J{x) = Jo{x), T{x) = Ji - ^Jo,a.(x) (120) 
In terms of these new variables, the second Hamiltonian structure takes the form 

{J(x),J(y)}2 = 2dJ{x-y) 
{T(x),J(y)}2 = J{x)dJ{x-y) 

{T{x),T{y)}2 = {T{x) + T{y))dJ{x-y) + Ul8{x-y) (121) 

Wc recognize this as the Virasoro-Kac-Moody algebra, which is the bosonic limit of the twisted 
N = 2 superconformal algebra. 
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Non-linear Schrodinger equation: 



Since the TB system reduces to the non-Unear Schrodinger equation, we can also find a Lax 
description for that system from the present one. We note that with the identification 



Jo - Qxloverq, Ji = qq 
the Lax operator for the TB system becomes, 

L = d-Jo + d-^Ji 

= d+^+d-'qq 
Q 

= q-^{d + qd-^q)q 
= GLG-^ 

G = q-^, L = d + qd-^q 



(122) 



where we have defined 



(123) 



(124) 



Namely, the two Lax operators L and L are related by a gauge transformation. The adjoint of this 
transformed Lax operator is determined to be 



C = V = -{d + qd-^q) 
It is straightforward to check that both the standard Lax equations 

di 



and 



at 

dt 



+ 



(125) 



(126) 



(127) 



give rise to the non-linear Schrodinger equation. (However, supersymmetry seems to prefer the 
second representation.) 



Furthermore, if we identify 



then, the standard Lax equations 



and 



give the mKdV equation. 



q = q = u 



dL 
dt 

dt 



(128) 
(129) 
(130) 



Thus, we see that the TB system is indeed a rich theory to study. 
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5 Supersymmetric equations: 



Given a supersymmetric integrable system, we can ask if there also exist supersymmetric integrable 
systems corresponding to it. It turns out to be a difficult question in the sense that the super- 
symmetrization turns out not to be unique and we do not yet fully understand how to classify all 
possible supersymmetrizations of such systems. Let me explain this with an example. 



Let us consider the KdV equation 



du 

'di 



Then, a supersymmetric generalization of this system that is also integrable is given by 

^ = 6uUx + Uxxx - ^Ipi^xx 

at 

dijj 

— = 3 (uV-)^ + V'xj;^ (131) 

Here represents the fermionic superpartner of the bosonic dynamical variable u of the KdV 
equation. It is easy to check that these equations remain invariant under the super symmetry 
transformations , 

Si/j = eu, Su = eipx (132) 

where e is a constant Grassmann parameter (fermionic parameter) of the transformation, satisfying 
e2 = 0. 



We can, of course, determine the supersymmetry charge (generator of supersymmetry) associ- 
ated with this system and it can be shown that the supersymmetry algebra satisfied by this charge 
is 

[Q,Q]+ = 2P (133) 

This can also be checked by taking two successive supersymmetry transformations in opposite 
order and adding them. There are several things to note from this system. First, the second 
Hamiltonian structure associated with this system, V2, is the superconformal algebra, which is the 
supersymmetrization of the Virasoro algebra. Second, just as this represents aN = 1 supersym- 
metric extension of the KdV equation, there also exists a second N = 1 supersymmetric extension. 



du 

'at 

— = Gui/jx + ipxxx (134) 

which is integrable. This second supersymmetric extension was originally discarded as being a "triv- 
ial" supersymmetrization, since the bosonic equations do not change in the presence of fermions. 

However, it generated a lot of interest after it was realized that it is this equation that arises in 
a study of superstring theory from the point of view of matrix models. Such a supersymmetric 
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extension now has the name -B supersymmetrization. Thus, we see that even at the level of N = 1 
super symmetrization, there is no unique extension of the integrable model. The problem becomes 
more and more severe as we go to higher supersymmetrizations. For = 2, it is known that there 
are at least four distinct "nontrivial" supersymmetrizations of the KdV equation that are inte- 
grable. Understanding how many distinct supersymmetrizations are possible for a given integrable 
equation, therefore, remains an open question. In addition, there are also fermionic extensions of 
a given integrable model (not necessarily super symmetric) that are also integrable and there does 
not exist any unified description of them yet. 



Lax description: 



Since the supersymmetric KdV equation (super KdV) is an integrable system, let us determine 
a Lax description for it. The simplest way to look for a Lax description is to work on a superspace, 
which is the natural manifold to study supersymmetric systems. 

Let us consider a simple superspace parameterized by {x,9), a single bosonic coordinate x and 
a single real fermionic coordinate 9, satisfying 6'^ = 0. Supersymmetry, then, can be shown to 
correspond to a translation, in this space, of the fermionic coordinate 6. On this space, one can 
define a covariant derivative 

which transforms covariantly under a supersymmetry transformation and can be seen to satisfy 

= d (136) 

On the superspace, a function is called a superfield and, since the fermionic coordinate is nilpotent, 
has a simple representation of the form (in the present case) 

<P{x,9) = i;{x) + eu{x) (137) 

The Grassmann parity of the components is completely determined by the parity of the superfield. 
For our discussion of the super KdV system, let us choose the superfield $ to be fermionic so that 
we can think of as the bosonic dynamical variable of the KdV equation and i/j as its fermionic 
superpartner. 

In terms of this superfield, the super KdV equations can be combined into one single equation 
of the form 

^ = (l>6$) + 3 (l>2 ($(£)$))) (138) 
It is now easy to check that if we choose, as Lax operator on this superspace, 

1 = 0^^ + 0^ (139) 



then, the Lax equation 

dL 

at 



(140) 
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gives the super KdV equation. The structure of this Lax operator and the Lax equation is, of 
course, such that they reduce to the KdV equation in the bosonic hmit, which is nice. 



It is worth making a few remarks about operators on the superspace. First, a pseudo-differential 
operator on this space is defined with powers of D. Correspondingly, various decompositions are 
done with respect to powers of D. Thus, we define 

super Residue P = sRes P = coefficient of 

super Trace P = sTr P = jdxdOsRes P (141) 

The conserved quantities, for the super KdV system, for example, are obtained as 

2n+l 

Hn = sTr L— (142) 

These are all bosonic conserved quantities and there is an infinite number of them. They all reduce 
to the infinite set of conserved charges of the bosonic integrable model in the bosonic limit. The 
Gel'fand-Dikii brackets can be generalized to this space as well and lead to the correct Hamiltonian 
structures of the theory. An interesting feature of the Lax description on a superspace is that 
the same integrable system can be described in terms of a Lax operator that is either bosonic or 
fermionic. Furthermore, in the supersymmetric models we can define non-local conserved charges 
from the Lax operator by taking, say, for example in the super KdV case, powers of quartic roots 
of the Lax operator. Finally, let me also point out that a supersymmetric -B system has fermionic 
Hamiltonians (conserved charges) with corresponding odd (fermionic) Hamiltonian structures. 



Super TB hierarchy: 



As we have seen, the TB hierarchy consists of two dynamical variables, Jo, Ji- Therefore, 
the supersymmetrization of this system will involve two fermionic partners, say ^Oj^i- Prom our 
experience with the super KdV system, let us combine them into two fermionic superfields of the 
forms 

$0 = V'o + OJo, «>i = V'l + OJi (143) 
With a little bit of algebra, we can check that if we choose the Lax operator 



then, we obtain 



(144) 
(145) 



which gives 

It is clear, therefore, that the Lax equation 

dL 
dt 



2{DHi)'^) + D-^ {iD^<^i) + 2D\'^i{D<^o))) (146) 



(147) 
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leads to consistent equations and gives 



dt 



{D^^i) + 2[D^{^i{D^o))) (148) 

The Lax operator as well as the Lax equation (and the equations following from it are easily seen 
to reduce to the TB equations in the bosonic limit. These, therefore, represent a supersymmetric 
extension of the TB hierarchy that is integrable. The higher order flows of the hierarchy are 
obtained from 

Once we have the Lax description of the system, we can immediately construct the conserved 
charges from 

Hn = sTi U' = jdxdesRes L", n = 1,2,3,- •• (150) 
Explicitly, we can construct the lower order conserved charges as 

= -Jdxde^i 
H2 = 2 j dxde{D^o)^i 

Hs = 3jdxd9(^{D^^o)-{D^i)-{D^of)^i (151) 

and so on. Thus, we see that the system has an infinite number conserved charges, which are in 
involution (follows from the Lax description) and, therefore, is integrable. 

Hamiltonian structures: 

Just like the bosonic system, the super TB equations also possess three Hamiltonian structures, 
although they are not necessarily local unlike the structures in the bosonic case. Let me only 
describe the first two here. 

Defining, as in the bosonic case, a generic Hamiltonian structure as. 



we note that 



as well as 



-2D-2D-^^iD-^ + D-\D'^<^o)D-^ - D{D^o) + D-^^iD \ 
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define the first two Hamiltonian structures of the super TB hierarchy. These structures have the 
necessary symmetry property and it can be checked using the method of prolongation that these 
structures satisfy the Jacobi identity. They give rise to super TB equation, say for example, through 



Although it is not obvious, the second Hamiltonian structure corresponds to the twisted N = 2 
superconformal algebra, which can be seen as follows. Let us look at the second Hamiltonian 
structure in the component variables. Let us define 



In terms of these variables, the nontrivial elements of the second Hamiltonian structure take the 
local form 



{Jo(.t), Jo(y)}2 = 2dJ(x-y) 

{Jo{x),Ji{y)}2 = da:{Jo6{x - y)) - dl6{x - y) 

{Ji{x),Ji{y)}2 = {Ji{x) + Ji{y))dJ{x - y) 

{Mx),i{y)}2 = ^Six-y) 

{Mx),m}2 = -mx-y) 

{Ji(x),e(y)}2 = (ax) + m)d.six - y) 

{Ji{x),ay)}2 = ax)dj{x~y) 

U»,e(y)}2 = -lji6{x-y) + ^dAJoS{x-y))-\dl6{x-y) (157) 



We recognize this to be the N = 2 superconformal algebra. 

N = 2 supersymmetry: 

Although the super TB system that we have constructed, naively appears to have A'' = 1 
supersymmetry, in fact, it does possess a N = 2 supersymmetry. This is already suggested by the 
fact that the second Hamiltonian structure of this system corresponds to the N = 2 superconformal 
algebra. Explicitly, this can be checked as follows. 

Let us note that the super TB equations, in terms of the redefined components, take the forms 




(155) 



v-l) 



(156) 




-KT — —Jo,xx + '^JqJq,x + 2Jl,a; 



X 



(158) 
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It is now straight forward to check that this system of equations is invariant under the following 
two sets of supersymmetric transformations, 



6Jo 
6Ji 

and 

6Jo = 2ef 
6Ji = 

k = (160) 



2eC 




(159) 



Nonlocal charges: 



As we have already noted, in supersymmetric integrable systems, in addition to the local bosonic 
conserved charges, we also have nonlocal conserved charges. Let us study this a little within the 
context of the super TB hierarchy. 



Let us recall that the Lax operator for the system is given by 

L = - {D^o) + D-'^^i 
and the infinite set of local, bosonic conserved charges are obtained as 

iJ„ = sTV L", n = 1,2,3,- •• 



(161) 



On the other hand, let us also note that in this system, we can also define a second infinite set of 
conserved charges as 

n = 1,2,3,- •• (162) 



2n-l 

Q 2n-i = sTr L 2 . 



There are several things to note from this definition. First of all, unlike the earlier set, these 
conserved charges are fermionic in nature. Second, they are nonlocal. Let me write down a few 
lower order charges of this set. 

Qi = - J dxd9{D-^^i) = - J dz{D-^^i) 

Qs = - j dz ^{D-^^if - $0^1 - {p-\{D^o)^^)) 

Q| = - j dz ^(D-i$i)=^- (5(D-2$i)$i_2$o$i-3(D$i)-(D-^$i)2)(D$o) 

+ (d-^ ({D^i)^i + ^liD^of - {D<^i){D^^o)))\ (163) 
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and so on. 



There are several things to be noted about these charges. First, as we have already mentioned 
and as can be explicitly seen, these charges are fermionic in nature and are conserved. Second, 
even though they are defined on the superspace, they are not invariant under super symmetry (it 
is the nonlocality that is responsible for this problem). This infinite set of charges satisfies an 
algebra which appears to have the structure of a Yangian algebra, which arises in the study of 
supersymmetric nonlinear sigma models, namely. 



{Qi,Qi}i 

2 2 


= 


{Qi , Qsji 

2 2 


= Hi 


{Qk,Q5.}i 

2 2 


= H2 


{Q3,Qs}i 

2 2 


= 2H2 


{QsjQbji 

2 2 


= '-Hs 
3 ^ 


2 2 


= 3H4 



7 

■ — J 
24 

5„ 



(164) 



and son on. The role of the fermionic nonlocal charges as well as the meaning of the Yangian 
algebra, however, are not fully understood. 



6 Dispersionless integrable systems: 



Let us consider again the KdV equation as an example. 

du g du d^'u 
dt dx dx^ 

As we have noted earlier, it is the linear term on the right hand side that is the source of dispersion 
for the solutions. Let us, therefore, get rid of the dispersive term, in which case, the equation 
becomes 

ft ^^ar.^ 

m = ^""Tx ^'''^ 

This is known as the Riemann equation and we sec that it corresponds to the dispersionless limit 
of the KdV equation. Given a nonlinear, bosonic equation, the systematic way in which the 
dispersionless limit is obtained is by scaling 

d d d d 

in the equation and then taking the limit e — 0. This leads to the dispersionless limit of the 
original system of equations. It is important to note here that the dynamical variable is not scaled 
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(although in supersymmetric systems, as we will see, it is necessary to scale the fermionic variables 
to maintain supersymmetry) . 



Let us recall that the Lax description for the KdV equation is obtained from the Lax operator 

L = d'^ + u 

through the Lax equation 



'di 



As we will now see, the Lax description for the dispersionlcss model is obtained in a much simpler 
fashion. Let us consider a Lax function, on the classical phase space, of the form 



L(p) =p'^ + u 



(167) 



Here p represents the classical momentum variable on the phase space and, therefore, this Lax 
function only consists of commuting quantities. However, we can think of this as consisting of a 
power series in p and formally calculate 



3 3 
p + -up 



(168) 



where the projections are defined with respect to the powers of p. It is now simple to check, with 
the standard canonical Poisson bracket relations. 



that the Lax equation 



{x,p} = l, {x,x} = 0{p,p} 
dL 



dt 



4{L, L2 I 



(169) 



leads to the Riemann equation, which is the dispersionless limit of the KdV equation. 



Thus, we see that given a Lax description of an integrable model in terms of Lax operators, the 
dispersionless limit is obtained from a simpler Lax description on the classical phase space through 
classical Poisson brackets. Let us also note that the conserved quantities of the dispersionless model 
can be obtained from this Lax function as (in this model) 

2n+l f 2n + l 

Hn = Ti L— = dx,Res L— , n = 0,l,2,--- (170) 

where "Res" is defined as the coefficient of the p^^ term. All of our discussions in connection with 
pseudo-differential operators carries through to this case where the Lax function has a polynomial 
structure in the momentum variables. (When I gave these talks, I had mentioned that the con- 
struction of the Hamiltonian structure from the Gel'fand-Dikii was an open question. Since then, 
this problem has been solved and we know now that these can be constructed rather easily from a 
Moyal-Lax representation of integrable models.) 
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Let me say here that dispersionless models encompasses a wide class of systems such as hydro- 
dynamic equations, polytropic gas dynamics, Chaplygin gas, Born-Infeld equation, Monge-Ampere 
equation, elastic medium equation etc, some of which show up in the study of string theory, mem- 
brane theory as well as in topological field theories. 



Let us study an example of such systems in some detail, namely, the polytropic gas dynamics. 
These are described by a set of two equations 

ut + uux + v'^~'^Vx = 0, 7 / 0> 1 

vt + {uv)x = (171) 

These equations are known to be Hamiltonian with 

= A (172) 

and 

so that we can write the polytropic gas equations as 

/ 5H \ 



Ut 
Vt 



V\ 1^ \ (174) 



5v 



In fact, this system has three distinct Hamiltonian structures, but I will not get into the details of 
this. 



Let us next consider a Lax function, on the classical phase space, of the form 

L = pT-Vn + Z"' \o P~^^~^^ (175) 
(7-1)^ 

Then, it is straight forward to check that the classical Lax equation 

gives rise to the equations for the polytropic gas. The higher order equations of the hierarchy are 
similarly obtained from 

f =.„((L-*)^^,.} (177) 

Thus, we see that the polytropic gas dynamics is obtained as a nonstandard Lax description on the 
classical phase space. (At the time of the school, this was the only Lax description for the system 
that was known. Subsequently, a standard Lax description has been obtained, which brings out 
some interesting connection between this system and the Lucas polynomials.) 
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Once we have the Lax description, we can, of course, obtain the conserved quantities from the 
"Trace". However, in this case, unHke in the case of Lax operators, we observe an interesting 
feature, namely, the residue can be obtained from expanding around p = or around p = oo. Thus, 
there are two series of conserved charges that we can construct for this system. Expanding around 
p = GO, we obtain 

Hn+i = Cn+iTr n = 0,1,2,- ■■ (178) 

where C„+i's are normahzation constants. Exphcitly, the first few of the charges have the forms 

Hi = J dxu 



Ho 



H. 



„7-l 



(7-l)(7-2)^ 



, / 1 o uv^-^ 

dx —U + —, r 

I 3! (7-i)(7-2)^ 



and so on. On the other hand, an expansion around p = leads to 

^„ = C'„TV L"+^, n = 0,1,2, • 
The first few charges of this set have the explicit forms, 



(179) 



(180) 



^0 

Hi 

H2 



-I 



dx V 
dxuv 

^"'2!"" + ^(^ 



(181) 



and so on. The two sets of conserved quantities can, in fact, be expressed in closed forms. Let us 
also note that if we define two functions as 



X = A T-i 



X 



A ^-1 



n + A 
2 

u + \ 



-,7-1 



(7 - I? 



2 u + X 
+ — - 



u + A 



1 

7-1 



1 

7-1 



(182) 



where A is an arbitrary constant parameter, it can be shown that these two functions generate the 
two sets of conserved quantities as the coefficients of distinct powers of A. Let me also note that 
the second Hamiltonian structure for the polytropic gas has the form 



+ ^,7-2^ du+{-f- 2)ud 
(7 — 2)du + ud dv + vd 



(183) 
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Dispersionless supersymmetric KdV: 



Let us recall that the super KdV equation can be described by the Lax operator 

L = D'^ + D^ 

and the Lax equation 



dt 



L 



Here, 

represents the covariant derivatiove on the superspace. In trying to obtain the dispersionless limit 
of this supersymmetric system, let us recall what we have learnt from the dispersionless limit of 
a bosonic model. We noted that the Lax operator goes over to the Lax function with d ^ p. 
However, our Lax operator, in the supersymmetric case, is described in terms of super covariant 
derivative D. Therefore, the natural question is what this object goes over to in the dispersionless 
limit. 

Let us note that the classical super phase space is parameterized by {x,6,p,p0) . Prom these 
we can define a variable 

U = -{pe + ep) (184) 
whose action on any phase space variable, through the Poisson brackets is 

{U,A} = {DA), {U,U} = -2p (185) 

Therefore, it would seem natural to let D — »• H in the dispaersionless limit. However, this leads to 
a serious problem. For example, we know that D"^ = d ^ p whereas H^ = since it is a classical 
fermionic variable. 

To analyze this problem a little more, let us recall that the dispersionless limit is obtained by 
scaling 

dt — > edti d ^ ed 
Therefore, since -D^ = d, consistency would require that we scale 

This implies that the fermionic coordinate needs to be scaled as 

e e-h (187) 
On the other hand, let us recall that the basic fermionic superfield of the theory is given as 

^{x,9)=i; + eu (188) 



34 



Since the dynamical variable u does not scale and scales, it follows that supersymmetry can be 
maintained under such a scaling only if 

'(p^e'^ip, $^e~5$ (189) 

This shows that, unlike in the bsoonic theory, in a supersymmetric theory, fermions scale in going 
to the dispersionless limit. With this scaling, we can go to the super KdV equation and note that, 
in the dispersionless limit, the equation becomes 

$t = 3D^ (190) 

which can be thought of as the super Riemann equation. 



Obtaining a Lax function and, therefore, a Lax description of a supersymmetric theory reamins 
an open question. However, through brute force construction it is known that the Lax function 

1 —2 —4 

L = p2 + + P_ ((D$)2 _ 2$$^) _ P—^Dmx (191) 

leads through the classical Lax equation (the projection is with respect to powers of p) 

§=4{L,(Li)^} (192) 

gives the dispersionless limit of the super KdV equation (super Riemann equation). It is worth 
emphasizing here that, although the Lax description for afcw supersymmetric dispersionless models 
have been constructed through brute force, a systematic understanding of them is still lacking. 



The conserved charges can be obtained from this Lax description immediately. Thus, 

Hn = C„TV = / ($(i5$)" - n$(£)$)"-^$^) (193) 

where n = 0, 1, 2, • • • and it is clear that these are bosonic conserved charges. The supersymmetric 
dispersionless model has the Hamiltonian structure 

P = -1 (3$D2 _^ (D$)D + 2(L>2$)^ (194) 

which can be recognized as the centerless superconformal algebra. With this Hamiltonian structure, 
it is straight forward to check that the conserved charges are in involution, 

{H^,Hrn} = ldz^-^V^-^=0 (195) 
which also follows from the Lax description of the system. 



Let me note in closing that this supersymmetric model has two infinite sets of nonlocal charges 
of the forms 

Gn = J dz^D-^'^)'' (196) 
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where n = 1,2,- ■■. (Note that Go = Hq.) These charges have been constructed by brute force, 
since it is not clear how to obtain nonlocal quantities from a classical Lax function. This remains 
an open question. Second, of the two sets of nonlocal charges, we see that F„ is fermionic while Gn 
is bosonic. Furthermore, all these conserved charges satisfy a very simple algebra, 

{Hn,Hm} = = {Fn, Hrn} = {Gn, Hrn} 
{Fn,Gm} = = {Gn,Gm} 

{Fn,Fm} = nmGn+m-2 (197) 



In connection with dispersionless supersymmetric integrable models, several questions remain. 
For example, it is not clear how to systematically construct the Lax description for them. It is not 
at all clear how nonlocal charges can be obtained from a classical Lax function. Neither is it clear 
what is the role played by these charges within the context of integrability. 

It is a pleasure to thank the organizers of the Swieca school for the warm hospitality. This work 
was supoported in part by US DOE Grant no. DE-FG-02-91ER40685. 
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